AD-A063  982 


UNCLASSIFIED 


WISCONSIN  UNIV-MADISON  MATHEMATICS  RESEARCH  CENTER  F/G  12/1 

EQUIVALENCE  OF  LINEAR  COMPLEMENTARITY  PROBLEMS  AND  LINEAR  PROGR— ETC(U) 
OCT  78  C W CRYER » M A DEMPSTER  DAAG29-75-C-0024 

MRC-TSR-1889  NL 


(H 

ADA063982 


technical  Nummary  Kep*r  1889 


EQUIVALENCE  OF  ^INEAR  COMPLEMENTARITY 
JROBLEMS  AND  yiNEAR  PROGRAMS  IN  ^ 
VECTOR  LATTICE  HILBERT  SPACES,  f I 


Mathematics  Research  Center 
University  of  Wisconsin-Madison 
610  Walnut  Street 
Madison,  Wisconsin  53706 

loctotaUBi  1978  / \ ] j 

(£[  J 

(Received  February  21,  1978)  "**1 

. 


Approved  for  public  rolooso 
Distribution  unlimited 


Sponsored  by 

U.  S.  Army  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


JJJ  300 

79  01 


National  Science  Foundation 
Washington,  D.  C. 

20550 


UNIVERSITY  OF  WISCONSIN  - MADISON 
MATHEMATICS  RESEARCH  CENTER 


r 


EQUIVALENCE  OF  LINEAR  COMPLEMENTARITY 
, PROBLEMS  AND  LINEAR  PROGRAMS  IN 

VECTOR  LATTICE  HILBERT  SPACES 

(1  2 3}  (4) 

. C.  W.  Cryer  ' ' and  M.  A.  H.  Dempster 

Technical  Summary  Report  #1889 
October  1978 
ABSTRACT 

* 

Let  X be  a vector  lattice  Hilbert  space  with  dual  X . Let  M be  a con- 

It  It 

tinuous  linear  mapping  of  X onto  X . Let  p,  q e X with  p > 0.  We  consider 
the  relationship  between  the  linear  complementarity  problem : Find  x e X such 
that  x > 0,  Mx  + q > 0,  <x,  Mx  + q)  = 0,  and  the  linear  programing  problem: 

/ Find  x c X which  minimizes  (x,p>  subject  to  x>0,  Mx+q>0. 
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SIGNIFICANCE  AND  EXPLANATION 


Many  free  boundary  problems  in  the  areas  of  fluid  mechanics,  porous  flow, 
elasticity,  and  plasticity  can  be  formulated  as  linear  complementarity  problems  in 
which  a differential  equation  (ordinary  or  partial)  must  be  solved  subject  to  the 
inequality  constraint  that  the  solution  be  non-negative?  roughly  speaking,  at  any 
point  the  solution  must  either  be  zero  or  satisfy  the  differential  equation.  We 
study  linear  complementarity  problems  which  can  be  reformulated  as  linear  programs 
in  which  a linear  functional  must  be  minimized  subject  to  inequality  constraints. 
The  reformulation  of  linear  complementarity  problems  as  linear  programs  offers 
two  advantages : 

(i)  It  suggests  alternative  numerical  methods  of  solving  the  problems. 

(ii)  For  the  problem  of  a cavitating  journal  bearing,  which  is  used  as  an  example, 
the  linear  program  requires  the  minimization  of  a linear  functional  which  is 
proportional  to  the  load  borne  by  the  bearing,  so  that  the  linear  program 
has  a physical  interpretation.  It  is  possible  that  the  linear  programs  for 
other  problems  will  also  have  physical  interpretations,  though  this  will 
have  to  be  determined  in  each  case. 
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EQUIVALENCE  OF  LINEAR  COMPLEMENTARITY  PROBLEMS 
AND  LINEAR  PROGRAMS  IN  VECTOR  LATTICE  HILBERT  SPACES 


C.  W.  Cryer*1'2'3*  and  M.  A.  H.  Dempster*4* 

1 . Introduction 

The  linear  complementarity  problem  in  real  n-dimensional  Euclidean  space  Rn  is:  Find 
n T 

x t R such  that  x > 0,  Mx  + q > 0,  and  x (Mx  - q)  “0,  where  M is  a given  real 

n x n - matrix  and  q is  a given  vector  in  Rn.  The  linear  programming  problem  in  Rn  is: 
n T 

Find  x e R which  minimizes  p x subject  to  x £ 0 and  Mx  + q > 0,  where  M is  a 

given  real  n x n matrix  and  p and  q are  given  vectors  in  Rn. 

Mangasarian  [1976]  showed  that,  under  certain  conditions,  each  solution  of  the  linear 

programing  problem  in  Rn  is  a solution  of  the  linear  conplementarity  problem  in  Rn. 

Mangasarian  [1977,  1977a]  has  subsequently  extended  this  work.  Related  work  is  due  to 
Cottle  and  Veinott  [1972],  More  [1971],  Tamir  [1973],  Cottle,  Golub,  and  Sacher  [1974], 
Cottle  and  Pang  [1976,  1976a],  Pang  [1977], 

Quite  independently,  and  often  not  very  explicitly,  the  relationship  between  certain 
infinite-dimensional  linear  programing  problems  and  linear  complementarity  problems  has 
been  noted  (Moreau  [1971] , Durand  [1968] , Lewy  and  Stampacchia  [1969] , Stampacchia  [1965] , 
Lions  and  Stampacchia  [1967]). 

Here,  we  consider  extensions  of  some  of  the  results  of  Mangasarian  to  infinite-dimen- 
sional spaces.  Apart  from  their  intrinsic  value,  our  results  provide  useful  ways  of  inter- 
preting, analysing,  and  solving,  linear  programing  problems  and  linear  complementarity 
problems  arising  in  physical  situations. 

The  paper  is  organized  as  follows.  In  section  2 we  dispose  of  some  preliminaries.  In 
section  3 we  introduce  the  linear  program  (LP),  the  dual  linear  program  (LD) , and  the  least 
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element  problem  (LE) . In  section  4 we  introduce  the  linear  complementarity  problem  (LC) , 
the  variational  inequality  (VI) , and  the  unilateral  minimization  problem  (or  quadratic  pro- 
gressing problem)  (UM) . In  section  S we  discuss  the  relationship  between  the  linear  program, 
the  least  eleswnt  problem,  and  the  variational  inequality.  In  section  6 we  discuss  in  detail 
a one-dimensional  problem,  finally,  in  section  7,  we  apply  our  results  to  the  problea  of 
lubrication  cavitation  in  journal  bearings. 


2.  Preliminaries 





X denotes  a real  Hilbert  space  with  norm  ||  • ||  and  dual  Y - X*.  The  evaluation  of 
a continuous  linear  functional  i e x*  at  a point  x t X is  denoted  by  <x,i>. 

It  is  assumed  that  X is  partially  ordered  by  a vector  ordering  >.  Let 

P - {x  t X : x > 0}  . 

Then  (Kelley  and  Namioka  [1976,  p.  224])  P is  a convex  cone  in  X with  vertex  at  the 

origini  that  is,  P P c p and  IP  c p for  all  non-negative  real  X.  We  assume  that  P 

is  closed,  x > y iff  x - y > 0,  that  is,  x - y e P. 

* * 

The  dual  cone  P c x is  defined  by 

• • * • 

P “ (x  c X : <x,x  > > 0 for  all  x e P)  . (2.1) 

* • * 

We  write  x >0  if  x e P . Since  P is  closed  it  follows  from  the  Hahn-Banach  theorem 

that  x » 0 iff  <x,x  ) > 0 for  all  x*  e P . 

It  is  also  assusied  that  X is  a vector  lattice  (Kelley  and  Namioka  [1976,  p.  229]). 
That  is,  for  all  x,y  e X,  there  exists  a unique  element  sup(x,y)  e X such  that 
sup(x.y)  > x and  sup(x,y)  > y;  furthermore,  if  z t X satisfies  z » x and  z > y then 
z > sup(x.y)  . The  assumption  that  X is  a vector  lattice  has  the  following  consequences. 
For  all  x,y  t X there  exists  a unique  element  infix, y)  such  that  x > infix, y)  and 
y > inf(x,y);  furthermore,  if  z e X satisfies  z < x and  z < y then  z < inf (x,y) . If 


x > y then 

sup(x,y)  - 

x,  and  if 

y > x then 

sup(x,y)  - yi  since 

sup(x,y)  is  unique,  it 

follows  that 

if  x > y 

and  y > x 

then  x • y. 

For  every  x « X, 

x “ sup(x,0)  - inf(x,0) 

so  that  X • 

P - p.  If 

0 - x + y 

where  x,y  e 

P then  x • y « 0; 

thus  0 is  an  extreme 

point  of  P, 

that  is. 

P is  a pointed  cone. 

* * * • 

M : X •*  Y ■ X denotes  a continuous  linear  transformation  with  adjoint  M i Y -*  X 

defined  by 


<x,  M*y*>  » <Mx,  y*>  . (2.2) 

Associated  with  M we  have  the  continuous  bilinear  operator  a : X x X -»  R1  defined  by 

a(v,u)  - (u,  Mv)  > (2.3) 

-3- 


(2.4) 


a is  symmetric  if  a(u,v)  ■ a(v,u),  and  coercive  if 

a(x.x)  > o||x||2  , 

for  some  real  strictly  positive  constant  a and  all  x t X. 


He  will  sometimes  impose  the  following  conditions  upon  a and  M: 


Condition  S.  If 

r « X 

and 

u,v  £ X are  such  that 

a(u,i|i) 

£ <i|i,r> 

and 

a(v,<|i)  > <<i,r>  for  all 

and  if 

w - inf (u,v) 

, then 

a(w,*) 

> <♦# r> 

for 

all  « P. 

Condition  Z.  If  u,v  £ P satisfy  inf(u,v)  - 0, 


then  a(u,v)  < 0. 


□ 


* * 

p and  q denote  elements  of  X . We  assume  throughout  that  p £ P . We  will  some- 
time assume  that  p is  strictly  positive,  that  is,  if  x £ P and  (x,p>  « 0 then  x « 0. 

Since 

a(u,<|*)  - <4>,r>  - <♦,  Mu  - r>  , (2.5) 

Condition  S may  be  rewritten  as  follows:  if  Mu  > r,  Mv  > r,  and  w * inf(u,v),  then 
Mw  > r.  If  -M  is  the  Laplacian  operator  V2  and  r - 0 then  Mu  > 0 means,  in  an  appro- 
priate sense,  that  -u  is  subharmonic.  In  this  case.  Condition  S reduces  to  the  well-known 
fact  that  the  infimum  of  two  superharmonic  functions  is  superharmonic.  There  is,  therefore, 
a close  connection  between  some  of  the  present  results  and  the  theory  of  subharmonic  func- 
tions (Rado  (1972],  Brelot  [1945,  1965],  Stampacchla  (1965],  Littman  [1963],  Moreau  [1971]). 

In  the  case  when  M is  a square  matrix,  Condition  Z is  equivalent  to  the  requirement 
that  the  off-diagonal  elements  of  M be  non-positive  - that  is,  that  M is  a z-matrix 
(Fiedler  and  Ptak  [1962]).  There  is,  therefore,  also  a close  connection  between  some  of  the 
present  results  and  the  theory  of  M-matrices  and  Z-matrices  (Plemnons  [1976]).  Condition  Z 
was  implicitly  used  by  Stampacchla  [1969,  p.  151]  with  the  conclusion  a(u,v)  <_  0 replaced 
by  a(u,v)  - 0. 

Conditions  S and  Z are  not  equivalent  because,  as  shown  in  section  2.1,  the  necessary 
and  sufficient  conditions  for  Conditions  S and  Z are  not  equivalent  in  the  case  of 
matrices . However,  we  do  have 
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Theorem  2.1 


Let  a be  coercive  and  satisfy  Condition  Z.  Then  a satisfies  Condition  S. 

Proof:  Let  u,v  £ X and  r e X satisfy  a(u,i|i)  > (<l>,r)  and  a(v,ip)  > <*,r>  for  all 

♦ > 0 . we  wish  to  show  that  if  w « inf(u,v)  then  a(w,i)>)  > <i|/,r>  for  all  i{i  > 0.  to  do 

so,  we  modify  an  argument  of  Stampacchia  [1965,  p.205] . 

Introduce  the  set  U c x which  consists  of  all  5 £ X satisfying  C > w.  U=p+w 
is  closed  and  convex.  From  the  fundamental  theorem  on  variational  inequalities  (Stampacchia 
(19641)  we  know  that  there  exists  n e U such  that 

a (n , z - n)  > (2  - n»  r>,  (2.6) 

for  all  z « U.  In  particular,  choosing  z - n + \|»,  we  see  that  a(q,^)  > (i|),r> 
for  all  il>  > 0.  The  theorem  will  therefore  be  proved  if  we  can  show  that  n « w. 

Set  £ » inf(r\,u)  £ U.  From  (2.6)  with  z » c, 

a(n,  C - n)  > (c  - n.  r>  . (2.7) 

On  the  other  hand  we  know  that  n - t > 0,  u - C > 0,  and  inf(n  - C,  u - C)  » 

inf(n,u)  -c»0.  Invoking  Condition  Z we  see  that 


a(c,  c - n)  * 

a (u. 

c - 

n) 

+ a (C  - u, 

C - n) 

- 

a (u, 

c - 

n) 

+ a (u  - C, 

n - C) 

< 

a(u,  C - 

n) 

t 

(C  - n» 

*> 

. 

Combining  (2.7)  and  (2.8)  we  find  that 

a(C  - n,  C - h)  < 0 . 

Since  a is  coercive  it  follows  that  c « inf(u,n)  “ n,  so  that  n $ u.  Similarly  n < v. 
Hence  n £ inf(u,v)  » w.  But  n e U so  that  n £ w.  We  conclude  that  n “ w. 

□ 

In  the  case  when  M is  a real  square  matrix,  it  is  readily  shown  from  Theorem  2.1.1 
below  that  if  a is  coercive  and  satisfies  Condition  S then  a satisfies  Condition  Z.  We 
do  not  know  whether  this  is  true  in  general. 

We  now  give  three  examples  of  spaces  and  operators  fitting  into  the  above  framework. 


2.1  Example  1 

* * n 

Let  X - Y » X • Y - R ; M - (m^)  *"  n 51  n reel  matrix;  and  p « (p^),  q - (q^)  » 

n-vectors.  Let  P be  the  set  of  vectors  in  Rn  with  non-neqative  components,  so  that  P 
* 

is  closed  and  P - P. 

P has  the  additional  important  property  that  it  has  non-empty  interior. 

Clearly,  p is  strictly  positive  iff  p^  > 0 for  all  i. 

It  is  readily  seen  that  Condition  Z is  satisfied  iff  w 5 0 for  i # j (that  is, 

M is  a Z-matrix) . 

Theorem  2.1.1;  Condition  S is  satisfied  iff  every  row  of  M has  at  most  one  strictly  posi- 
tive coefficient,  that  is,  MT  is  pre-Leontief  (Cottle  and  Veinott  11972,  p.  244])  . 

T 

Proof;  We  first  observe  that  M is  pre-Leontief  iff  each  row  k of  the  inequality 
Mu  ^ r can  be  written  in  the  form 


il 

c u > r.  + 7 d.u. 

S S - k ^11 


where  d.  > 0 for  all  j;  c > 0;  and  where  the  dependence  upon  k of  c and  d.  has 
J - s « s j 

been  suppressed. 

T 

First  let  us  assusie  that  M is  pre-Leontief  and  that  Mu  > r,  Mv  > r.  Then  inequality 
(*)  holds  for  u,  and  a similar  inequality  holds  for  v.  Since  d^  > 0 we  have  that  if 
w - inf (u, v)  then 


c w - inf(c  u , c v ) , 

s s s s s s 


n 

> r ♦ l d.  inf (u  ,v  ) , 

j-1  ■'  3 J 

n 

■ rk  ♦ l djwj  * 
j-i  3 3 


so  that  Condition  S is  satisfied. 


Now  let  us  assusw  that  Condition  S holds  but  that  M is  not  pre-Leontief.  Then  there 
is  a row  k of  N with  at  least  two  positive  coefficients,  m^  and  m^  say.  Thus,  the 
k-th  row  of  the  Inequality  Mu  > r takes  the  form 


« 


\sus  i rk  - "ktut  - ^ -kj-j  ' 

«nd  a similar  inaquality  holds  for  v.  Set  ug  » 1/m^,  ufc  - vg  - 2ug,  vt  - 2ut< 

and  Uj  “ v^  ■ 0 otherwise.  Finally,  set  w « inf(u,v)  and 

rj  ■ mini X -vt ' X "jtv  ' 

for  all  j.  Then  Mu  > r,  and  Mv  > r.  But, 

n n 

lx  "Vi  ■ "ksus  + "ktvt  + "kjuj  ' 

- r - 1 i*a,t 

rk  1 ' 

so  that  the  inequality  Mw  > r does  not  hold. 


In  conclusion,  we  note  that  if  the  problems  in  Examples  2 and  3 below  are  discretized, 
the  resulting  finite  difference  or  finite  element  matrices  usually  satisfy  Conditions  S and 


2 . 2 Example  2 


Let  X ■ Hg(fl),  where  ft  is  a bounded  domain  (open  connected  set)  in  Rn,  and  H^(ft) 

is  the  Sobolev  space  of  once-differentiable  functions  vanishing  on  3ft  (Adams  11975)). 

* -1 

■Then  Y ■ X ” H (ft).  Let  M be  the  linear  self-adjoint  operator. 


(Mu)  (t)  - -l  jf-O^ft)  ^r-(t))  , t e ft  , 


(2.2.1) 


with  coefficients  a^(t)  which  are  continuously  differentiable,  where  the  indices  i and 
j are  simned  from  1 to  n.  It  is  assumed  that  -M  is  uniformly  elliptic,  so  that 


Iaij(t)ei«j  > aU|2  , t £ ft 


(2.2.2) 


for  all  C - (ti>  £ R , and  some  constant  a > 0. 

Every  x £ H*  has  a representation  as  a measurable  function  x(t),  and  any  two  such 

representations  of  x differ  only  on  a set  of  measure  zero.  He  write  x > 0 if  x(t)  > 0 

a. a.  (almost  everywhere).  P - {x  « X : x > o)  is  clearly  convex. 

-7- 


I 


To  show  that  P is  closed,  let  {x  } be  a sequence  of  points  in  P which  converges 

n 

to  x e H*.  Then  x (t)  converges  to  x(t)  in  l?  (SI)  , from  which  it  follows  that 
0 n 

x (t)  -*■  x(t)  a.e.  Hence,  x(t)  > 0 a.e.  so  that  x € P. 

n ” 

hJ  is  a vector  lattice:  if  x,y  £ then  the  functions 


sup(x,y)  (t)  « sup(.-(t),  y(t)),  , t e fJ  , 


inf(x,y)(t)  “ inf(x(t),  y(t)),  , t £ (2  , 


(2.2.3) 


are  representations  of  elements  in  HQ.  (Lewy  and  Stampacchia  [1969,  p.  169]  prove  that 
H1(0)  is  a vector  lattice,  and  their  proof  can  be  readily  adapted  to  the  present  case.) 

Another  very  useful  property  of  H*  is  that  if  x £ hJ  and  F is  a measurable  subset 
of  0 on  which  x(t)  is  constant  then  (Lewy  and  Stampacchia  [1969,  p.  169]), 


/ | grad  x(t) | dt  = 0 . 

F 


(2.2.4) 


As  defined  in  (2.2.1),  the  operator  M can  only  be  applied  to  functions  u which  are 
twice  differentiable.  Let  a:H*  * H*  ■*  R 1 be  the  symmetric  coercive  bilinear  operator  de- 
fined by 


0 i j 


(2.2.5) 


We  extend  the  domain  of  definition  of  M by  regarding  M as  the  mapping  from  X = HQ  to 

* -i 

its  dual  space  X « H defined  by 


<v, Mu  > - a (u,v) , for  all  u,v  £ H 


(2.2.6) 


The  standard  theory  of  elliptic  operators  allows  us  to  assert  that  M is  uniquely  defined 

1 * -1 

by  (2.2.6)  and  that  M is  a homeomorphism  of  X ■ HQ  onto  X “ H (Lions  and  Magenes 
[1972,  p.207]).  In  particular,  M is  an  open  mapping. 

Theorem  2.2.1 

M satisfies  Conditions  S and  Z. 

Proof:  To  prove  Condition  Z,  let  u,v  £ P and  inf(u,v)  » 0.  Let  u vanish 
on  P c 0 and  v vanish  on  G c 0.  Then,  using  (2.2.4),  we  conclude  that 


a(u,v)  « If  aij(t)  , 

n i j 


n ' ' * ' 1 

Wij 3t7  3t 


*/•«<  «£-&■«♦  /^(U^^at-o 


j 


G-F 


so  that  Condition  Z is  satisfied. 

Stampacchia  [1965,  p.  205]  proves  that  Condition  S is  satisfied. 

□ 

* 

We  conclude  by  making  some  additional  remarks  about  inequalities  in  X and  X . 

(a)  Let  x e Hg(ft).  Then  x > 0 in  the  sense  of  H^fQ)  if  there  exists  a sequence 

{W  } of  functions  f e C^(£ 5)  which  satisfy  <p  (t)  > 0 in  {l  and  which  converge 
in  m m = 

to  x in  H1(!l)  (Lewy  and  Stampacchia  [1969,  p.  155]).  If  x > 0 in  the  sense 
of  H1(0)  then  it  follows  immediately  that  x(t)  > 0 a.e.  . Conversely,  let  x « 
satisfy  x(t)  > 0 a.e.  . If  x denotes  the  extension  of  x to  Rn  obtained  by  set- 
ting x(t)  « 0 for  t i 0,  we  know  that  x t H^(Rn)  (Adams  [1975,  p.  57]) . The 
averaged  functions  xh  are  smooth  and  non-negative,  and  they  converge  to  x in 

(Rn)  (Adams  [1975,  p.  52]).  If  ^ - x.  1(1  then  >p.  ■*  x in  0,  and  we  can 

n n n 

conclude  that  x > 0 in  the  sense  of  H1(n) . We  have  thus  shown  that  if 
x e Hg(fl)  then  x > 0 in  the  sense  of  H^fO)  iff  x(t)  > 0 a.e.  . This  is  of  im- 
portance to  us  because  Stampacchia  and  his  colleagues  use  > 0 in  the  sense  of  H^(fi). 

(b)  If  r e H ^ is  non-negative  then  there  exists  a non-negative  Borel  measure  u 
such  that 

(W,r)  - f W dli 
ft 

for  every  * e c”(fl)  (Schwartz  [1973,  p.  29]).  This  is  a very  elegant  charac- 
terization of  non-negative  functionals;  unfortunately,  it  is  difficult  to  apply 
because  its  use  involves  measure  theory. 

2.3  Example  3 

This  example  is  the  special  case  of  Example  2 when  0 - (a,b)  c R1 . 
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M.1  the  properties  of  Example  2 remain  valid.  There  is  one  additional  property  which 
is  very  important:  Every  x e H^(a,b)  can  be  represented  as  a function  x(t)  which  is 
absolutely  continuous  on  [a,b]  and  vanishes  at  tne  endpoints.  (Smirnov  [1964]). 


a,  x 


1 £ x ( HQ(a,b)  then  x « x^  + x^  + xf  where  x^  vanishes  outside  a neighborhood  of 

vanishes  outside  a neighborhood  of  b,  and  x has  compact  support  in  (a,b) . The 

1 in 


averaged  functions  x 


(h) 


of  x are  smooth,  have  compact  support,  and  converge  to  x , 


as  h -*•  0.  We  can  construct  similar  approximations  x|h*  and  x 


.(h) 


to  x.  and  x by 
l r 


first  translating  and  then  averaging.  If  x > 0 a.e.  then  f.  » x. 

■ n l 

verges  to  x and  we  see  that  x > 0 in  the  sense  of  H^(a,b) . 


(h) 


♦ x(h)  ♦ x<h> 
m r 
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3.  The  linear  program,  the  dual  linear  program,  and  the  least  element  probli 


With  the  notation  of  section  2,  the  linear  program  (LP)  is: 

(LP)  Minimise  <x,p>  subject  to  Mk  + q > 0 

x<P 

The  dual  program  (LOP)  which  is  (formally)  dual  to  LP  is: 


(LOP) 


where 


Maximise  (-q,y  ) subject  to  -My  + p > 0 
y *P 


• • * 

{y  « y : y >0}  , 


- * * * * * •- 

* (y  € y : <u,y  > > 0 for  all  u € P } . 

If  x is  a solution  of  LP  and  y is  a solution  of  LET  then, 

• * * * 

(x.p)  - <-q,y  ) - (x,  - M y + p>  ♦ (Mx  + q,  y ) > 0 , 


(3.1) 


(3.2) 


(3.3) 


(3.4) 


so  that  the  value  of  LP  is  always  greater  than  or  equal  to  the  value  of  LDP.  In  particular, 

• • * 

if  <x,p)  * (q.y  ) • 0 for  some  feasible  x and  y then  x and  y are  optimal.  It  may, 

however,  occur  that  the  two  values  are  never  equal,  in  which  case  there  is  a duality  gap. 

Since  x is  reflexive  we  know  (Dunford  and  Schwarts  (1966,  p.  66])  that  there  is  an 

*•  # 

isometric  isomorphism  k which  maps  X onto  X ■ y and  which  is  defined  by 


Let 


where  y t X (not  Y ) , so  that 


(x,x  ) - <x  ,oc) 

y*  - *y  . 

(-q,y*>  - <-y,q> 


(3.5) 

(3.6) 

(3.7) 


Me  assert  that  y >0  iff  y t P.  First  asauaw  that  y c P.  Then,  for  any  u e P 

* * • 

<y,u)  > 0,  so  that  y ^ 0.  On  the  other  hand,  suppose  that  y > 0 but  that  y f P.  Then, 

since  the  singleton  (y)  is  compact  and  the  cone  P • {x  t X : x > 0)  is  closed  and  convex, 
these  two  sets  can  be  separated  (Dunford  and  Schwarts  (1966,  p.  417]).  That  is,  there 
exists  a linear  functional  f t Y and  constants  c > 0 and  c such  that 
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« 


Using  the  properties  of  P we  conclude  that  c « 0,  so  that  ftp  and  <y,f)  < But 


<y,f>  « <f,y  >>  0,  and  we  have  a contradiction 
Finally,  for  any  u e X , 


(definition  of  M ) 


(equation  (3.6)) 


(definition  of  k) 


where  we  define  the  linear  operator  M : X -*  X « Y by 


Thus 


Summing  up,  we  see  that  y satisfies  LDF  iff  y » icy  where  y solves 


and  we  will  take  this  to  be  the  dual  of  LP  in  our  further  work 


Since  X is  partially  ordered,  we  nay  also  consider  the  least  element  problem  (LE) 


Find  x e P such  that  Mx  + q > 0 and  x < u for  every  u £ P satisfying  Mu  + q > 0 


were  two  solutions  we  would  have  x,  < x 


LE  has  at  most  one  solution,  for  if  x.  and  x, 


and  x_  < x.  which  implies  that  x. 


there  exists  a very  satisfactory  theory  for  LP  and  LD 


and  Mangasarian  11976]  used  this  as  the  starting  point  for  his  study  of  the  relationship  be 
tween  LP  and  LCP  (the  linear  complementarity  problem) . LE  has  also  been  studied  in  the 


finite  dinensional  case  (Cottle  and  Veinott  (1972] ) 


The  case  when  X is  infinite  dinensional  is  nuch  more  difficult.  It  is  usually  assumed 


[1974,  p.  66],  that  the  Arrow-Hurwicz  constraint  qualifica- 


tion is  satisfied,  namely  that  there  exists  u £ P such  that  Mu  + q is  an  interior  point 


■ 


I 


4.  The  linear  complementarity  problem,  the  variational  inequality,  and  the  unilateral 
minimi ration  problem. 

The  linear  complementarity  problem  (LC)  is  as  follows:  Find  x e P such  that 

(LC)  Mx  + q > 0,  <X,Mx  + q)  - 0 . (4.1) 

“ 

The  variational  inequality  (VI)  is:  Find  x « P such  that 
(VI)  a(x,  v - x)  + <v  - x,  q)  > 0 , (4.2) 

for  all  v £ P. 

If  a is  synnetric  then  the  unilateral  minimisation  problem  (or  quadratic  programing 
problem)  (UM)  is:  Find  x e P such  that 

(UM)  J(x)  < J(u),  for  all  u e P (4.3) 

where 

J(u)  - a(u,u)  + 2<q,u)  . (4.4) 

The  basic  result  on  variational  inequalities  is  due  to  Stampacchia  11964] : if  a is 
coercive  then  there  exists  a unique  solution  to  VI. 

The  connection  between  VI  and  UM  was  also  observed  by  Stanpacchia  [1964] : if  a is 
symetric  and  coercive,  then  VI  is  equivalent  to  UM. 

The  relationship  between  VI  and  LC  was  noted  independently  by  a number  of  workers  in- 
cluding Lions  and  Stampacchia  11969,  p.  172],  Karamardian  11971],  More  11971].  The  basic 
result,  which  we  prove  for  the  convenience  of  the  reader  is: 

Theorem  4.1 

LC  is  equivalent  to  VI. 

Proof : Assume  first  that  x solves  LC.  Then  for  any  v t P 

a(x,  v - x)  + <v  - x,  q>  - <v  - x,  Mx  ♦ q)  , 

“ (v,  Mx  + q>  > 0 , 

so  that  x solves  VI. 

Now  assume  that  x solves  VI.  Then  setting  first  v - 0 and  then  v - 2x  we  see 
that  (— x,  Mx  ♦ q>  > 0 and  <x,  Mx  + q) £ 0,  from  which  we  conclude  that  <x,  »t  ♦ q)  • 0. 
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. the  relationship  between  the  linear  program,  the  least  element  problem 


[ 


complementarity  problem. 

Theorem  S . 1 

If  a is  coercive  and  satisfies  Condition  Z,  then  LE  has  a solution,  namely  the  unique 
solution  of  VI. 

Proof;  The  proof  is  a modification  of  proof  of  Stampacchia  [1969,  p.  151]  who  implicitly 

used  Condition  Z in  the  special  form:  if  u,v  « P and  inf(u,v)  » 0 then  a(u,v)  = 0. 

Let  u be  the  unique  solution  of  VI  so  that  u e P and 

a(u,  v-u)  + (v-u,  a > >0 

for  all  v r P. 

In  particular,  choosing  v ■ u + w for  any  w e p we  conclude  that  Mu  + q > 0. 

Now  let  w be  any  element  such  that  w e P and  Mw  + q > 0.  We  assert  that  w > u. 

TO  see  this,  let  C “ min(u,w)  e X,  so  that  w - C > 0 and  u - C > 0.  Furthermore, 

inf(w  * {,  u • {)  ■ inf(w,u)  - t » 0. 

Then 

a(u  - C,  u - O - [a(C,  C - u)  + <t  - u,  q>]  - [a(u,  5 - u)  + - u,  q>]  , 

< la(C,  C - u)  + <t  - u,  q)]  , 

m 

because  u satisfies  VI.  But 

a(£»  C - u)  + (C  - u,  q)  » a(w  - C,  u - C)  + (a(w,  C - u)  + <C  - u,  q)]  , 

< 0 , 


because  the  first  term  on  the  right  is  nonpositive  by  Condition  Z and  the  second  term  is 

nonpositive  since  Mw  + q > 0 and  5 - u < 0. 

Combining  the  above  inequalities  we  see  that  a(u  - C,  u ~ O £ Remembering  that 

a is  coercive  we  conclude  that  u “ C.  Thus,  w ^ C ■ u so  that  u is  a solution  of  LE. 

□ 

Theorem  5.2 

(i)  If  x solves  LE  then  x solves  LP . 

(ii)  If  a satisfies  Condition  S and  p is  strictly  positive,  then  LP  has  at  most 
one  solution. 

-16- 


I 


i 


(iii)  If  • satisfies  Condition  S,  p is  strictly  positive,  and  x solves  If  then 
x solves  LE. 

Proof:  (i)  is  obvious.  To  prove  tii),  let  xx  and  x?  be  two  solutions  of  LP.  By  Condi- 
tion S,  c - inf(x1.x2)  e P satisfies  Me  + q > 0 and  <C,p>  * (x^p).  since  x is  optimal, 

<C'P>  " (*X.P>  and  we  conclude  that  C - xy  Similarly,  5 - x2»  so  that  Xj  - 

To  prove  (iii),  let  u « P satisfy  Mu  - q > 0.  Set  C - inffu.x).  Then  M<  + q > 0 

and  <C,p>  “ <x,p)  so  that  C “ x.  Hence,  u > x and  x solves  LE. 

□ 

Remembering  that  if  a is  coercive  and  a satisfies  condition  Z then  a satisfies 
Condition  S (Theorem  2.1)  we  find  that 
Theorem  S.3 

If  a is  coercive  and  satisfies  Condition  Z,  and  if  p is  strictly  positive,  then  LP, 
LE,  VI,  and  1C  all  have  the  same  unique  solution. 

□ 

Theorem  5.4 

Assume  that  x solves  VI,  that  y solves  LD,  that  (x,p)  + <y,q)  - 0,  that  a is 
symmetric  and  coercive  and  satisfies  Condition  Z,  and  that  p + q > 0. 

Then  y > x. 

Proof:  Set  w«inf(x,y).  Then 

a(x  - «r,  x - v)  - a(x  - y,  x - w)  + a(y  - w,  x - w) , 

< a(x  - y,  x - v)  , 

since  y - w > 0,  x - w > 0,  and  inf (y  - w,  x - w)  - 0.  But, 
a(x  - y,  x - w)  - a(x,  x - w)  - a (y,  x - w)  , 

- a(x,  x - w)  - a(x  - w,  y)  , 

- a(x,  x - w)  + <x  - w,  - My)  , 

- [a(x,  x - w)  + <x  - w,  q>]  - <x  - w,  p ♦ q>  ♦ 

+ <x  - w,  - My  + p>  . 
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The  first  term  on  ths  right  Is  negative  because  x solves  VI.  The  second  term  is  negative 
because  p ♦ q t P and  x - w c P.  The  third  term  is  zero  because  the  equality 

0 - <x,p>  ♦ <y.q>  » <x,  - My  + p>  + <y,  Mx  + q> 

Implies  that  (x,  - Sy  + p>  - 0 and  hence,  since  0 < w < x,  that  <w,  - My  + p>  - 0. 

Combining  the  above,  we  conclude  that  a(x  - w,  x - w)  < 0 so  that  x ■ w.  Then 
y >»■  inf (x,y)  - x. 

□ 

It  may  be  observed  that  if  x solves  LP,  y solves  LD,  (x,p)  + <y,q)  - 0,  and 
y > x,  then  we  have  that 

0 < <x,  Hx  ♦ q>  < <y,  Mx  + q)  - 0 ; 

that  is,  x solves  1C. 
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6.  A one-dimena.onal  problem 


Ha  consider  a special  case  of  Example  3 (section  2.3)  i X-H*(0,2), 

2 

Min(x,p)  ■ / lx(t)dt  subject  to  x(t)  > 0 a.e., 

0 

and  Mx  + q ; -x(t)  + (t  - 1)  ^ 0 , 

with  the  corresponding  dual  problem 

2 

*“<y*_q)  “ - / (t  - l)y(t)dt  subject  to  y(t)  > 0 a.e., 

0 " 

-My  + p = y(t)+l>0  . 

The  inequality  -X  + (t  - 1)  > 0 is  interpreted  in  the  sense  that 

2 

<*,  NX  ♦ q>  - /(i(t)e(t)  + (t  - l)*(t))dt  > 0 , 

0 


(6.1) 


(6.2) 


(6.3) 


All  non-negative  ( hJ(0,2),  and  the  inequality  y + 1 > 0 is  interpreted  in  the 

This  problem  was  chosen  because  it  is  a simple  problem  with  the  same  general  structure 
••  the  problem  for  a cavitating  journal  bearing  which  is  discussed  in  the  next  section. 

There  is  a straightforward  procedure  for  obtaining  possible  solutions  of  such  one-dimen- 
sional problems)  these  solutions  can  then  be  verified  a-postiori.  Me  assume  that  x(t)  > 0 for 
0 < t < t and  x (t)  • 0 for  T < t < 2,  where  t is  an  unknown  constant  corresponding 
to  the  free  boundary  (the  point  t - t) . If  x also  satisfies  LC  then  <x,  -X  ♦ (t  - 1) > -0, 
so  that  -3i(t)  ♦ (t  - 1)  ■ 0 for  0 < t < t.  The  general  solution  of  the  equation 
-Si  ♦ (t  - 1)  ■ 0 is 

x(t)  - A ♦ Bt  ♦ i (t  - l)3  . (6.4) 


Using  the  conditions  x(0)  - x(t)  - 0 to  determine  the  constants  A and  B we  find 


x (t)  - t(t  - x) [-3  ♦ t ♦ t]/6 


(6.5) 
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To  determine  x we  note  that  the  condition  -x  + (t  - 1)  >0  implies  that  for  all 
smooth  non-negative  f « (0, 2) , 


(W,  Mx  + q>  » / [x  i*  (t  - l)*]dt  , 
0 


t 2 

» / [xp  + (t  - l)W)dt  + / (t  - l)t*>  dt  * 
0 x 

1 X 2 

<f  r + / i-*>  + (t  - D^]  dt  + / <t  - dv>  dt 
J°  0 X 


x(x-)v(x)  + / (t  - 1)*  dt  , 


This  is  only  possible  if  x > 1 (so  that  t - 1 > 0 for  t < lx, 2])  and  x(x-)  > 0.  But, 
x(t)  > 0 for  t < x and  x(x)  - 0 so  x(x-)  < 0.  We  conclude  that  x(x-)  » x(x+)  - x(x) 
• 0.  The  condition  x(x)  - 0 leads  to  an  algebraic  equation  for  x,  namely. 


i(X)  - x(-3  + 2XJ/6  - 0 ; 


thus,  x “ 3/2  and 


xlt)  - tit  - 3/21/6  , 0 < t < 3/2  , 


0 , 3/2jt<2 


is  our  trial  solution. 

Using  (6.6)  and  (6.7)  we  see  that  x is  such  that  x > 0,  rx  + (t  - 1)  > Q,  and 
(x,  Wb  ♦ q)  * 0,  so  that  x is  a solution  of  LC.  InvoRing  Theorems  4.1  and  5.3,  we  con- 
clude that  x is  the  unique  solution  of  LP. 

We  now  consider  the  determination  of  y.  Since  (x,  y + 1)  « 0,  it  follows  that 
y (t)  ♦ 1 - 0 when  x(t)  > 0,  that  is,  when  0 < t < x.  On  the  other  hand,  since 
(y,  -5  ♦ (t  - 1))  - 0,  it  follows  that  y (t)  » 0 when  -x  + It  - 1)  >0,  that  is,  when 
x « t < 2.  We  conclude  that  y(t)  +1-0  for  0 < t < 3/2  and  y(t)  - 0 for  3/2  < t < 2. 
Solving  this  boundary  value  problem  we  obtain 
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y (t)  - tl-2t  ♦ 3]/*>  0 < t < 3/2  , 

•0  , 3/2  < t < 2 . 


(6.8) 


The  condition  y > 0 is  seen  to  be  satisfied. 
Direct  compytation  yields 


2,  9 2 

<x,p>  - / x (t) dt  - r-rr-  ■ - / (t  - l)y(t)dt  - -<y,q>  . 

0 128  0 


(6.9) 


The  solutions  x(t)  and  y(t)  are  plotted  in  Figure  6.1.  We  note  that  y ^ t as 
proved  in  Theorem  5.4. 


Figure  6.1:  x(t)  and  y(t) 

It  is  possible  to  give  two  justifications  for  the  free  boundary  condition  x(T) 
2 

Firstly,  if  x t Hq(0,1),  as  is  often  the  case,  then  $(t)  is  continuous  so  that 
4(t)  « Hi*)  • 0.  Secondly,  a reasonable  interpretation  of  the  condition 
-x(t)  ♦ (t  - 1)  > 0 is  that 
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7 . Lubrication  cavitation  of  journal  bearings 


A large  number  of  physical  problems  can  be  formulated  as  linear  complementarity  problems 
in  which  a differential  equation  (ordinary  or  partial)  must  be  solved  subject  to  the  inequal- 
ity constraint  that  the  solution  be  non-negative;  roughly  speaking,  at  any  point  the  solution 
must  either  be  zero  or  satisfy  the  differential  equation  (Cryer  [1977,  in  preparation], 

Duvaut  and  Lions  11972] ) . The  reformulation  of  such  linear  complementarity  problems  as  linear 
programs  has  two  advantages:  (i)  it  suggests  alternative  methods  of  solving  the  problems; 
and  (ii)  it  sometimes  provides  a physically  meaningful  interpretation.  As  an  example  of  such 
linear  complementarity  problems  we  consider  here  the  problem  of  a cavitating  journal  bearing. 

A journal  bearing  consists  of  a circular  cylinder  (the  journal)  which  is  rotating  inside 
a support  structure  (the  bearing) . The  narrow  gap  between  the  journal  and  the  bearing  is 
filled  with  a thin  film  of  lubricating  fluid.  Various  geometries  are  possible.  In  Figure 
7.1  we  show  a partial  journal  bearing  of  finite  length.  The  term  'partial'  refers  to  the 
fact  that  the  journal  is  not  completely  enclosed  within  the  bearing,  and  is  partially  ex- 
posed to  the  atmosphere. 
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It  is  required  to  determine  the  pressure  x of  the  lubricant,  and  the  load  W borne 
by  the  bearing.  Because  the  gap  between  the  journal  and  the  bearing  is  very  narrow,  the 
siaplifications  of  lubrication  theory  can  be  applied.  In  particular,  it  is  assumed  that  the 
pressure  does  not  vary  across  the  gap,  so  that  the  problem  becomes  a two-dimensional  problem 
in  the  rectangular  dastain  Q ■ ABCDEF  in  the  8 z -plane  (Figure  7.2). 


Figure  7.2:  The  domain  fl 


The  lubricant  flows  in  from  a reservoir  along  the  entry  edge  AT  and  flows  out 
through  the  ends  ABC  and  DEF  as  well  as  through  the  exit  edge  CD.  At  all  these  points 
the  lubricant  is  in  contact  with  the  atmosphere,  and  if  the  pressure  is  normalized  so  that 
atmospheric  pressure  is  zero,  then  the  boundary  conditions  are  that  x - 0 on  30.  That  is. 
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* « x - H*(n)  . 


(7.1) 


The  lubricant  occurs  In  both  liquid  and  gaseous  phases.  It  is  assumed  that  the  lubri- 
cant vaporizes  when  the  pressure  is  zero,  so  that  the  inequality  x > 0 must  be  satisfied 
everywhere . If  the  pressure  is  greater  than  zero  then  the  lubricant  is  in  the  liquid  phase 
and  satisfies  the  simplified  form  of  the  Navier-stokes  equations  known  as  Reynolds'  equation. 
After  introducing  dimensionless  variables,  the  equation  takes  the  form  (Pinkus  and  Sternlicht 
(1961)): 


MX  + q S 


_L  ,h3  3x 

36  <h  36  ) 


+ a 


_3_ 

3z 


<h‘ 


— ) 
3z  1 


+ 0 

de 


(7.2) 


where  a is  a positive  constant,  and  where  h • h(8)  is  a given  function  which  is  propor- 
tional to  the  width  of  the  gap. 

On  the  free  boundary  T,  the  interface  between  the  liquid  and  gaseous  phases,  the 
boundary  conditions  are 

x “ 0,  3x/3n  « 0,  on  T , (7.3) 

where  3/3n  denotes  the  normal  derivative. 

In  the  engineering  literature  (Pinkus  and  Sternlicht  [1961] ) the  problem  is  formulated 
mathematically  as  a classical  free  boundary  problem:  Find  x and  r such  that  x satis- 
fies (7.2)  subject  to  the  boundary  conditions  (7.1),  and  (7.3).  However,  in  a large  number 
of  papers  in  the  engineering  literature,  beginning  with  the  work  of  Christopherson  (1941) , 
numerical  approximations  have  been  obtained  in  a completely  different  way:  equation  (7.2) 
is  replaced  by  finite  differences,  and  the  resulting  system  of  algebraic  equations  is  solved 
ss  a finite-dimensional  linear  complementarity  problem  (Cryer  (1971))  which  may  be  considered 
as  a discretization  of  the  infinite-dimensional  linear  complementarity  problem 

x > 0,  »c  + q > 0,  (X,  Mx  + q)  - 0 . (7.4) 

Ne  may  thus  take  (7.4)  as  the  starting  point  for  a mathematical  analysis  of  the  problem. 

The  problem  is  a special  case  of  Example  2 (section  2.2),  and  it  follows  from  Theorem  S.3 
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that  there  exists  a unique  solution  x < H*(0)  of  LP,  LE,  VI,  and  LC. 

In  the  engineerinq  literature,  there  has  been  jcas  discussion  of  an  appropriate  varia- 
tional principle  for  the  problem  (Christopherson  [1957] ) . The  formulation  as  a variational 
inequality  leads  to  two  useful  variational  principles: 

(1)  Since  a is  symmetric,  the  problem  is  equivalent  to  the  unilateral  minimisation 
problem 

Inf  J(v)  - a(v,v)  + 2<v,q>  . 

v>0 

m 

(2)  Por  any  strictly  positive  function  p(6,s),  the  problem  is  equivalent  to  the 
linear  progressing  problem 

Min  <x,p>  = / x(6,s)p{e,s)de  ds  , 

a 

subject  to  x > 0,  Mx  + q > 0.  In  particular,  if  -n/2  * ®p  * ®D  * w/2  (see  Figures 
7.1  and  7.2),  then  p - cos  8 > 0 and  <x,p)  is  the  load  w borne  by  the  bearing  in 
^ vertical  direction  (Figure  7.1).  That  is,  the  solution  x minimises  the  vertical 
load. 
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